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This  memorandum  describes  a  method  of  calculating  the  equivalent 
circuits  of  discontinuities  in  coaxial  lines  caused  by  abrupt  changes 
of  radius  in  one  or  both  of  the  coaxial  conductors.  The  theory  is 
fully  worked  out  for  discontinuities  of  the  inner  conductor  alone; 
the  implementation  of  the  final  equations  on  a  large  computer  is 
described;  and  the  performance  of  the  resulting  program  is  discussed, 
including  the  effects  of  truncation  error  and  extrapolation  to 
infinity.  Finally,  the  use  of  the  theory  in  connection  with  coaxial¬ 
line  reflection  and  transmission  standards  is  described,  and  a  simple 
error  analysis  of  these  standards  is  presented,  which  shows  that  (at 
present)  their  performance  is  calculable  to  an  accuracy  limited  only  by 
the  precision  to  which  they  can  be  made. 
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INTRODUCTION 


This  memorandum  describes  a  practical  method  of  calculating  the  equivalent 
circuits  of  coaxial-line  discontinuities  of  the  types  shown  in  Figs  la-le.  The 
method  used  can  be  applied  to  many  similar  problems  in  one-conductor  or  two- 
conductor  waveguide  systems,  as  described  by  Harrington!  1] ;  but  Harrington's 
treatment  is  completely  general  and  gives  no  feeling  for  the  computational 
details  which  are  essential  when  the  method  is  applied  to  a  particular  case. 
Accordingly,  the  treatment  presented  here  discusses  fully  two  particular  cases 
which  are  important  in  practice,  and  uses  them  to  illustrate  the  general  theory, 
thus  complementing  Harrington's  exposition. 

The  method  applies  in  principle  to  any  problem  involving  semi-infinite 
"generalised  cylindrical  waveguides"  joined  at  a  plane  at  right-angles  to  their 
common  axis.  (A  generalised  cylindrical  waveguide  is  defined  by  translating  a 
simple  closed  plane  curve  at  right-angles  to  its  own  plane,  as  shown  in  Figs  2.) 
The  behaviour  of  the  electromagnetic  field  in  the  neighbourhood  of  the  dis¬ 
continuity  plane  in  such  a  system  can  be  characterised  by  an  equivalent  circuit 
located  at  that  plane;  the  discussion  which  follows  will  show  this,  establish 
that  the  form  of  the  equivalent  circuit  is  a  single  shunt  element  (in  general) 
and  obtain  an  expression  for  its  value.  The  numerical  aspects  of  calculating 
the  value  will  be  described  also,  as  will  the  practical  importance  of  the  work 
(for  standards  of  coaxial-line  impedance)  and  the  factors  affecting  accuracy. 

2  PRELIMINARIES  ON  FIELDS  AND  MODES 

The  first  stage  in  the  treatment  of  any  particular  case  is  to  express  the 
fields  in  the  guides  in  terms  of  the  permitted  modes  for  that  case.  The  theory 
of  modes  in  generalised  cylinders  has  been  fully  discussed  elsewhere  (eg  [  lj  , 

[2]  and  [ 3] ) ,  so  this  section  will  only  present  a  summary  of  the  important  facts. 

If  a  region  of  space  (such  as  the  interior  of  a  generalised  cylinder)  is 
linear,  isotropic,  homogeneous  and  non-conducting,  and  contains  no  free  charges 
or  current-carrying  elements,  any  electromagnetic  field  within  it  may  be  written 
as  the  sum  of  two  partial  fields,  one  of  which  has  its  E_  vector  transverse  to 
an  arbitrary  fixed  direction  in  space  (a  transverse  electric  or  TE  field)  while 
the  other  has  its  H  vector  transverse  to  the  same  arbitrary  fixed  direction  (a 
transverse  magnetic  or  TM  field);  and  each  partial  field  can  be  derived  from  a 
single  scalar  function  which  satisfies  the  wave  equation.  When  we  deal  with 
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generalised  cylinders  having  a  common  axis  it  is  convenient  to  make  that  axis 
the  arbitrary  fixed  direction  and  call  it  the  z-direction;  then  fields  which 
vary  in  time  as  exp(jwt)  can  be  derived  from  scalar  functions  ¥  satisfying 
Helmholtz*  equation 


where  k  is  the  free-medium  angular  wavenumber  corresponding  to  o>  and  V^,  is  the 
two-dimensional  Laplacian  operator  transverse  to  z.  These  functions  V  satisfy 
boundary  conditions  which  are  independent  of  z,  and  so  eqn  (1)  is  separable, 
giving 

(V*  +  k'2)<Jj  =  0  (2a) 


with 


H !  =  H 


d2C  2 

72  '  * 

dz 


(2b) 

(2c) 


and 


2 

Y 


(2d) 


The  boundary  conditions  on  ^  become  boundary  conditions  on  ip;  and  it  turns  out 

that  eqn  (2a)  admits  non-trivial  solutions  only  for  particular  non-negative  real 

2  2  . 
values  of  k*  ,  and  corresponding  values  of  y  from  eqn  (2d),  each  with  its  own 

2  2 

corresponding  field.  Each  of  these  sets  (k*  ,  y  and  corresponding  field) 
constitutes  a  mode,  and  any  possible  field  within  a  generalised  cylinder  can  be 
expressed  as  a  linear  combination  of  all  the  modes  for  that  cylinder.  (This 
does  not  conflict  with  the  division  into  TE  and  TM  fields,  because  each  mode 
can  be  classified  as  TE  or  TM. ) 


Power  is  freely  propagated  only  by  the  fields  of  those  modes  which  admit 

oscillatory,  rather  than  exponentially-decaying,  solutions  to  eqn  (2c),  that  is, 

2 

those  which  have  negative  values  of  y  .  There  are  only  a  finite  number  of  these, 
and  the  number  depends  on  the  frequency;  in  one-conductor  systems  at  sufficiently 
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low  frequencies  there  are  none  at  all.  However,  in  two-conductor  systems  there 

2 

is  always  one  mode  having  k*  =  0,  which  will  propagate  power  at  any  frequency. 
This  mode  belongs  to  the  TM  class,  but  has  special  properties  and  is  often 
considered  to  be  in  a  class  of  its  own;  because  its  j£  and  H  vectors  are  both 
transverse  to  the  z-axis,  it  is  called  transverse  electromagnetic  (TEM) . 

3  ROTATIONALLY-SYMMETRIC  SYSTEMS 

It  is  possible  to  write  down  the  modes  which  can  exist  in  systems  like 
those  in  Fig  1,  from  the  discussions  in  [1]  and  [2],  Alternatively,  we  can 
derive  them  directly  by  using  rotational  symmetry  about  the  z-axis  instead  of 
the  translational  symmetry  along  it. 

Subject  to  the  restrictions  introduced  already.  Maxwell's  equations  for  a 
field  varying  as  exp(jwt)  take  the  form  (see  [2]) 


curl  E  +  jwpH  “  0 

(3a) 

div  H  *  0 

(3b) 

curl  H  -  jmeE  =  0 

(3c) 

div  _E  =  0 

(3d) 

in  the  usual  notation. 

We  introduce  cylindrical  polar  co-ordinates,  expand  the  vector  operators 
in  eqns  (3)  and  enforce  rotational  synmetry  by  setting  to  zero  all  the  partial 
derivatives  with  respect  to  </>,  giving  (see  Fig  3) 


3E 

"  aT*  +  juijHp  =  0 

(4a) 

H  (pE;  +  =  0 

(4b) 

3H  3H 

_P  jueE^  -  o 

(4c) 

,  .  3H 

~  ■—  (PH  )  +  —  -  0 

p  dp  P  3z 

(4d) 

-  -r - jweE  ■  0 

3z  J  p 

(5a) 

7  h  (PHJ  “  -  0 

P  <f>  2 


(5b) 


(5c) 


9E 
_ P 

9z 


9E 


w  +  juyH<! 


o 


1  9  9E2 

p  Ip  «V  *  v 


(5d) 


On  substituting  eqns  (4a)  and  (4b)  into  eqns 
(4d)  is  satisfied  identically  and  that  (4c)  takes 


(4c)  and  (4d),  we  find  that 
the  form 


(PE  . )  -  i  (PE  ) 
<P  P  op  <P 


(pE.)  +  k2(PE  ) 


0 


(6) 


with 

k2  ■  a )2pe  (6a) 

as  usual.  A  similar  treatment  of  eqns  (5)  leads  to 
2  2 

^-9  <PH*)  ~  7  |r  (PH.)  +  2—  (pH  )  +  k2(pH  )  =  0  (7) 

3p2  ♦  p  9p  ♦  3Z2  *  ♦ 

The  boundary  conditions  to  be  used  with  eqns  (6)  and  (7)  follow  from  the 
requirements  that  the  tangential  electric  field  on  any  perfectly-conducting 
surface  must  vanish.  Using  eqns  (5a)  and  (5b),  we  deduce 

<»V  ■  »  <8a>  i;  <»V  ‘  0  (8b) 

where  n  denotes  a  normal  to  the  boundary  surface  (or  surfaces).  Eqns  (4),  (6) 

and  (8a)  specify  a  TE  field,  being  obtained  by  solving  (6)  subject  to  (8a) 

and  the  dependent  components  H  and  H  being  given  by  (4a)  and  (4b);  similarly 

p  z 

eqns  (5),  (7)  and  (8b)  specify  a  TM  field  having  components  ,  Ep  and 
4  CYLINDRICAL  SYMMETRY 

To  apply  these  results  to  the  systems  represented  in  Figs  1,  we  begin  by 
solving  eqn  (7)  using  eqns  (2).  Writing 

(pH  )  =  iKpK(z)  (9a) 

$ 
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►  A’  • 


we  obtain 


-  y2£  »  0  and  y 2  *  k'^  -  k^,  as  in  (2),  and 

dz 


I  ii 

P  dp 


k'2i> 


0 


(9b) 


2  2 
If  k'  *  0,  eqn  (9b)  has  the  solution  ip  -  Ap  +  B,  where  A  and  B  are  arbitrary 

constants;  otherwise  the  solution  is 


=  (k'p)(A'J1(k’p)  +  B'Y^k'p)) 


(9c) 


where  A'  and  B'  are  arbitrary  constants  and  and  Y^  are  Bessel's  functions 
of  the  first  and  second  kinds  of  order  1,  as  may  be  shown  by  substitution. 

Eqn  (6)  gives  similar  solutions  for  (pE^) ,  but  we  shall  see  that  these  are  not 
actually  needed. 

The  systems  of  Figs  lare  made  up  of  right  circular  cylinders  butted 

together.  Accordingly,  everywhere  except  at  the  plane  of  discontinuity  the 

£ 

boundary  conditions  (eqns  (8))  take  the  form  (pE^)  *  0  and  (pH^)  =  0  at  one 
or  two  particular  values  of  p,  depending  on  whether  the  system  contains  one  or 
two  conductors;  and  if  there  is  only  one  conductor  the  field  components  must  be 
finite  at  p  =  0. 


It  is  usual  to  operate  coaxial  systems  under  conditions  such  that  the  TEM 
mode  is  the  only  one  which  can  freely  propagate  power;  from  the  boundary  con¬ 
ditions  and  the  special  solution  we  find  that  this  mode  is  given  by  (pH^)  ■  ipQ  =  Be 

(no  non-trivial  special  solution  for  (pE,)  is  possible,  and  the  solution  for 

$ 

(pH^)  also  vanishes  if  there  is  only  one  conductor).  The  TEM  mode  is  a  TM  field, 
and  systems  like  those  of  Figs  1  (with  rotationally-symnetrical  discontinuities 
as  well  as  conductors)  have  independent  TM  and  TE  fields.  It  follows  that  only 
a  TM  field  can  exist  in  the  systems  of  Figs  1  when  they  are  worked  so  that  the 
TEM  mode  is  the  only  one  freely  propagating;  any  TE  field  travelling  inwards 
from  infinity  would  decay  to  zero  before  encountering  the  plane  of  discontinuity. 
Consequently,  if  we  obtain  the  permitted  modes  for  the  systems  of  Figs  1  by 
treating  them  as  generalised  cylinders,  we  must  reject  all  TE  modes,  as  well  as 
all  TM  modes  which  are  not  rotationally-symmetrical. 
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If  these  rejected  modes  are  considered,  it  turns  out  that  at  least  one 
of  them  begins  to  propagate  power  freely  at  a  frequency  at  which  all  the 
retained  modes  are  still  cut  off.  (See  [2]  for  a  full  discussion.)  Because 
of  this,  there  are  two  "critical  frequencies"  for  any  system  of  the  type  con¬ 
sidered  here;  a  lower  critical  frequency  above  which  at  least  one  (TE)  mode 
propagates  freely,  and  an  upper  critical  frequency  above  which  at  least  one  of 
the  retained  modes  propagates  freely.  (Conveniently,  the  TM  mode  with  the 
lowest  cut-off  frequency  is  one  of  the  retained  modes,  otherwise  -  since  the 
TE  modes  and  the  rotationally-asymmetrical  modes  are  rejected  for  different 
reasons  -  there  would  be  three  critical  frequencies.)  The  analysis  which 
follows  is  formally  valid  below  the  upper  critical  frequency,  but  will  only 
hold  between  the  two  critical  frequencies  provided  that  all  TE  modes  have  been 
filtered  out  before  one  of  the  discontinuities  represented  in  Figs  1  is 
encountered.  Expressions  for  the  critical  frequencies,  after  [2],  will  be  given 
later. 


5  DEVELOPMENT  OF  THE  PROPER  MODAL  EXPANSION 


We  are  now  ready  to  construct  the  modal  expansion  of  the  allowed  field  in 
any  of  the  systems  of  Figs  1;  the  one  shown  in  Fig  lb  will  be  used,  because 
its  practical  importance  was  the  reason  for  undertaking  the  work  described  here. 
In  this  case  the  boundary  condition 


0 


dip 

dp 


must  be  satisfied  at  p  =  and  p  =  R  on  the  A-side  of  the  discontinuity  plane, 
and  at  p  =  r  and  p  =  R  on  the  B-side.  The  TEM  mode  already  satisfies  this 
condition;  for  the  TM  modes,  using  eqn  (9c)  and  the  relations 


^  (xJ x (x) )  =  xJq(x)  and  ^  (xY^x)) 


xY  (x) 
o 


» 


we  obtain 

(k'p)(A'J  (k’p)  +  B'Y  (k’p))  =  0  (10) 

o  o 

to  be  satisfied  at  the  proper  values  of  p. 

On  the  A-side,  the  only  non-trivial  values  of  k'  satisfying  eqn  (10)  are 
those  which  satisfy 
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0 


(10a) 


J  (k'r.)Y  (k’R)  -  J  (k’R)Y  (k'r.)  = 

O  A  O  O  O  A 

From  Watson  [4],  the  roots  of  eqn  (10a)  are  real,  positive  and  distinct;  they 
will  be  denoted  by  k^,  where  i  is  a  positive  integer  and  kA1  <  kA2  <  kA3  ...  . 
Similarly,  on  the  B-side  we  have  a  set  of  kfi^  (kg^  <kB2<  kg^  •••)  satisfying 

Jo(k'rB)Yo(k’R)  -  JQ(k'R)Yo(k’rB)  =  0  (10b) 


The  boundary  condition  also  determines,  for  each  allowed  k',  the  proportions 
of  the  combination  of  J^(k'p)  and  Y^(k'p)  in  the  corresponding  ip.  We  must  take 
care  that  this  proportionality  does  not  become  infinite,  so  we  write  the 
allowed  tp  functions  on  the  A-side  as 


"Ai 


Ai  Ai 


'( 


A;,-(k*,p)|J1(kAiP)  + 


(ID 


with  a  similar  equation  on  the  B-side;  the  degenerate  form  of  this  is  then 
AAi^kAiP^Jl^kAip^  which  is  the  proper  function  for  the  degenerate  form  of 
Fig  lb  (namely  Fig  la).  The  boundary  condition  dtj^/dp  =  0  at  p  =  rA  and  p  =  R 
then  gives,  taking  Aj^  +  0, 


(k  .  r. 
Ai  A 

(k.  . r. 
Ai  A 


) 
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J  (k..R) 
o  Ai 

Y  (k..R) 
o  Ai 


(12) 


We  see  that  eqn  (10a)  is  the  consistency  condition  associated  with  eqn  (12). 
A  similar  result  holds  good  on  the  B-side. 

We  now  define  the  following  six  families  of  functions: 

Z  (k. . p)  =  J  (k  . p)  +  (B7A')..Y  (k.  . p) 
o  Ai  o  Ai  Ai  o  Ax 


Zl<kAiC>  ■  Jl<kAio)  *  <B’'A'>AiW’> 


(13a) 


22<kAi°>  -  J2(kA.p)  .  (B'/A')AiY2(kA.p) 


VkBip)  ’  VkBip)  *  (B'/A'Jg.Y^OcgjP) 


zl<kBlP)  '  -W>>  *  (B’/A'^t^.p) 


(13b) 


Z2(kBiP)  -  J2(kB.p)  .  (B'/A’)B.Y2(kB.p) 
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noting  that  these  Z  functions  must  always  have  arguments  of  the  form  (k. .p)  or 

Al 

(k„.p)  so  that  the  (B'/A')  factor  is  defined.  From  eqns  (11),  (12)  and  (13) 

D1 

we  obtain 


Z  (k..r.)  -  Z  (k..R)  *  Z  (k_.r.)  *  Z  (k„.P)  =  0 

o  Ai  A  o  Ai  o  Bi  B  o  Bi 


(14) 


for  all  i, conveniently  expressing  the  boundary  conditions,  and 


^Ai 


A;i(kAip)zi(kAr> 


(15a) 


with  a  similar  equation  on  the  B-side. 

To  complete  the  construction  of  tne  mode  functions,  we  refer  to  eqn  (9a) 
and  the  equations  immediately  following  it;  these  give 


(pfyAi  =  *Ai(p)<Ai(;s) 


(15b) 


with  £^(z)  satisfying 


d\ i  2  . 

2  *  YAiCAi  "  0 

dz 


and 


Y.U 


2  2 
k  -  k 
*Ai  *A 


(15c) 


From  eqn  (6a)  we  have 


,,2  2 

kA  =  “  UAEA 


(15d) 


in  an  obvious  notation;  and,  as  always,  similar  equations  hold  on  the  B-side. 

The  functions  c(z)  for  the  modes  which  do  not  propagate  power  freely  must  give 

fields  which  decay  as  we  move  away  from  the  discontinuity  plane;  so,  taking  z 

to  be  zero  at  that  plane  and  positive  on  the  B-side,  we  have  c  functions  of  the 
YAiz  ~YBiz 

form  e  and  c  (taking  y^,  Yg^  all  positive).  No  arbitrary  multiplying 

constant  is  needed  with  these  functions,  since  there  is  already  one  in  each 
and  only  the  products  are  physically  significant.  For  the  TEM  mode  we  must 
include  fields  freely  propagating  in  both  directions  and  on  both  sides  of  the 


discontinuity  plane,  and  from  eqn  (15c)  we  have  =  -  jk^  or  ±  jk-g*  The 


complete  field  is  a  linear  combination  of  all  these  possible  mode  fields  (by 
[  1]  and  [2]),  so,  if  we  write 


A! . 
Ai 


jooe, 


kA.y..  i 
Ai  Ai 


^i 


jue 


8. 


kBiYBi  1 


£ 


(16) 


in  terms  of  the  material  constants  y  ,  c^,  yB  and  Eg,  then  eqns  (15)  give 


H,  (z  <0) 


-  —  a  ( 

v  °V 


-jk.z  jk  z' 


A"  _  J"A  .  . 

e  -  PAe  ) 


Zl(kAi°) 


H  (z  >  0) 


J-e  (- 
V  o  V 


JkB2  r  “jkBZ\  . 
e  +  rBe  ) 


P®  ■ 


-YB:Z 

e  Z. (k_  .p) 


'1  Bj 


(17a) 


The  signs  and  multiplying  factors  here  have  been  chosen  for  convenience 

later.  These  equations  hold  for  the  degenerate  case,  Fig  la,  as  well  as  in  the 
case  of  Fig  lb,  except  that  F  vanishes  when  the  inner  conductor  is  truncated. 
To  complete  the  representation,  we  use  eqns  (5)  and  (17a)  to  obtain 


E  (z  <  0) 
P 


■  i«  ( 

p  o  \ 


jV  .  r  ->kAZ 
e  +  T.e 

A 


p) 


i-1 


E  (z  >  0) 
P 


-  S  (< 
P  °  \ 


JkBZ  r  -jkBZ 

e  +  rBe 


p)  • 


(17b) 


3-1 


s 


and 


E2  <*<»)  -  2(-^)°i'Y  ',Zzo<kAi‘» 

00 

Ez  (z  >  0)  -  B.z  TBjZ  Zo(kB.») 


To  obtain  eqns  (17c)  we  use,  as  before,  the  relations 
(xJ1(x))  =  xJq(x)  and  ^  (xY1(x))  =  xYq(x) 

together  with  eqns  (13)  to  deduce  that 


(17c) 


dOTTo y  ((kAip)Zl(kAip))  =  (kAip)Zo(kAip) 
Ai 


with  a  similar  equation  in  kfij  .  Eqns  (14)  and  (17c)  confirm  that  our  represen¬ 
tation  satisfies  all  the  boundary  conditions  introduced  so  far. 

rA  and  can  be  regarded  as  voltage  reflection  coefficients  if  Ep,  or 
some  quantity  proportional  to  it,  is  taken  as  the  measure  of  "voltage".  We 
will  refer  to  this  later. 


6  GENERALISED  MODES  AND  CRITICAL  FREQUENCIES 

The  most  general  set  of  modes  in  coaxial  systems  includes  others  not 
derived  here,  depending  on  Bessel  functions  of  higher  orders  and  exhibiting 
rotational  asymmetry.  These  are  of  no  significance  in  the  theory  of  the  equi¬ 
valent  circuit  of  the  discontinuity,  but  they  must  be  taken  into  account  in 
determining  the  critical  frequencies  referred  to  previously. 

The  upper  critical  frequency  is  derived  from  eqns  (10) ,  since  it  turns  out 

that  the  TM  mode  with  the  lowest  frequency  of  free  propagation  is  rotational ly- 

symmetrical  (it  is  the  one  usually  denoted  TMQ^) ;  it  begins  to  propagate  freely 

at  the  lower  of  the  two  frequencies  v  ,  vD  given  by 

A  D 

2-aV^  "  kA  =  kAl  (18a) 


11 


2  IT 


VB  V 


¥b 


=  kT 


B1 


(18b) 


using  the  notation  already  introduced.  The  lower  critical  frequency,  however, 
relates  to  a  rotational ly-asynmetrical  mode,  the  so-called  TE^  mode  (which  is 
actually  a  pair  of  degenerate  modes  distinguished  only  by  their  planes  of 
polarisation);  and  to  characterise  this  mode-pair  we  must  replace  eqns  (10)  by 

Jl(k'rA)Yi(k'R)  '  Ji(k'R>*i(fc,rA)  =  0  (19a) 


and 


J|(k,rB)Y^(k’R)  -  J^(k'R)Y’(k’rB)  =  0  (19b) 

where  the  prime  on  the  Bessel  functions  denotes  differentiation.  The  lowest 
values  of  k*  satisfying  these  equations  are  then  used  in  eqns  (18)  to  find  the 
critical  frequency. 


Good  approximate  solutions  of  these  equations  may  be  deduced  from  the 
discussion  in  [2]  and  are  as  follows: 


lower  critical  frequency 


upper  critical  frequency 


smaller  of  1/(tt  Vv*  (R  +  rA}) 
and  1 / ( n  (R  +  rB)) 

smaller  of  1/(2 ^AeA  (R  “  rA)) 


and  1/(2  (R  “  rB)) 


If  we  take  R  =  3.5  mn,  r.  =  rn  =  1.52  ram,  u  *  Up  *  u  »  e.  *  =  e  (the 

parameters  of  a  7  nun-bore  precision  air-spaced  coaxial  line) ,  these  approximate 
formulae  give  the  critical  frequencies  as  19.0  and  75.7  GHz,  the  exact  values 
being  19.4  and  75.1  GHz  to  the  same  number  of  significant  figures.  It  is  worth 
noting  that  changing  the  radius  of  the  inner  conductor  causes  the  critical 
frequencies  to  move  in  opposite  directions,  so,  when  uA  ■  Vg  and  ■  eB  (as 
is  usually  the  case)  the  lower  critical  frequency  is  fixed  by  r and  the  upper 
one  by  rB>  (r  is  greater  than  rfi;  see  Fig  lb.)  These  conclusions  need 
modifications  in  detail  for  the  systems  of  Figs  Ic  and  Id,  but  the  principles 
are  the  same. 
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For  the  system  of  Fig  la,  eqns  (10b)  and  (19b)  reduce  to  the  much  simpler 

forms 


JQ(k'R)  =  0  (20a)  and  Jj(k’R)  =  0  (20b) 

The  lowest  values  of  k'  satisfying  these  equations  are  very  nearly  (2.405/R) 
and  (1.841/R)  respectively,  from  which  we  obtain  approximate  critical  frequencies 
as  follows: 

lower  critical  frequency  =  smaller  of  1  / ( it  (R  +  r)) 

and  0.920/(tt^pb£b  R) 

upper  critical  frequency  =  smaller  of  1/ (2  (R  -  r)) 

and  1.202/(7T^pBeB  R) 

If  we  take  the  same  dimensions  and  material  constants  as  before,  the  approxi¬ 
mate  critical  frequencies  are  19.0  and  32.8  GHz,  the  exact  values  being  19.4 
and  32.8  GHz  to  the  same  number  of  significant  figures.  The  interesting  point 
here  is  the  large  reduction  in  the  upper  critical  frequency;  this  makes  the 
behaviour  of  the  equivalent  circuit  change  with  frequency  considerably  more 
rapidly  than  it  does  in  the  ordinary  coaxial  system  (Fig  lb) . 

7  THE  FIELDS  IN  THE  PLANE  OF  DISCONTINUITY 

We  are  now  ready  to  begin  the  derivation  proper.  In  the  introduction  it 
was  promised  that  this  would  be  presented  in  a  way  which  illustrated  the  proper 
approach  to  other  similar  problems;  so  we  note  that  equations  corresponding  to 
eqns  (17)  will  exist  for  any  problem  involving  generalised  cylinders  and  the 
first  stage  of  any  derivation  is  always  to  find  the  corresponding  equations 
representing  the  fields.  (Various  important  special  cases  are  discussed  in 
[2J.)  Next,  we  determine  the  transverse  fields  at  the  discontinuity  plane,  and 
apply  the  boundary  condition  that  they  must  be  continuous  (since  the  tangential 
components  of  _E  and  H  are  continuous  across  any  surface  not  supporting  a  surface 
current).  There  is  one  other  boundary  condition  to  be  applied:  the  transverse 
components  of  E  must  vanish  on  those  portions  of  the  discontinuity  plane  which 
are  parts  of  the  perfectly-conducting  cylinders  (rB  <  p  <  rA  in  the  system  of 
Fig  lb,  0  <  p  <  r  in  that  of  Fig  la,  for  example).  All  this  is  readily 
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generalised;  in  our  particular  case,  writing  and  for  the  values  of  Ep 
and  H  at  the  discontinuity  plane  z  *  0,  we  obtain 


£  =  -  a  (1  +  r.)  +  7  a.Z, (k. .p) 

p  p  o  A  i  1  At 


i=l 


7  6  (l  +  r  )  + 

P  O  D 


Zyi 

j=i 


(kBjp) 


(rA  <  P  <  R) 


(rg  <  p  <  R) 


(rB  <  P  <  rA)  ) 


(21) 


V.  -  —  a  (1  -  r.) 
'VA  n.P  O  A 
A 


S(-3*)« 


(kAip) 


(rA  <  p  <  R)  ^ 


6 


nBP  o 


(rfi  <  p  <  R) 


K  (22) 


These  equations  completely  express  the  boundary  conditions  at  z  *=  0,  since 
in  the  present  problem  Ep  and  are  the  only  transverse  components  of  the  E 
and  H  fields. 


8  CURRENT  AND  VOLTAGE  IN  THE  PRINCIPAL  MODE 


The  next  stage  in  the  development  is  to  obtain  expressions  for  the  "current" 
and  "voltage"  at  the  transverse  plane  z  =  0.  In  the  present  case,  since  we  are 
dealing  with  a  two-conductor  system,  natural  meanings  can  be  attached  to  these 
quantities  because  the  transverse  field  distribution  arising  from  a  TEM  mode  is 
essentially  independent  of  frequency.  If  we  make  use  of  this  fact  to  reduce  the 
problem  to  an  electrostatic  one,  eqn  (3a)  becomes 

curl  jJ  =  0  which  implies  that  I!  =  -  grad  V.  From  this  we  obtain 

dp  (23a)  taking  the  inner  conductor  as  positive 
with  respect  to  the  outer  one,  where 
[  1  means  "the  TEM  part  of". 


outer 

radius 


voltage  = 


./ 


inner 

radius 
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In  electrostatic  terms,  this  corresponds  to  making  the  inner  conductor  the  line 
and  the  outer  one  the  return.  The  current  in  the  inner  is  then  positive,  and 
is  related  to  the  magnetic  field  it  produces  by  the  most  general  electrostatic 
form  of  eqn  (3c) ,  which  is  (see  [  2] ) 

curl  fl  =  J.  Applying  Stokes'  theorem  (see  [2])  to  this,  taking  $  as  a 
simple  closed  curve  around  the  inner  conductor  lying  in 
the  plane  z  =  0  and  taking  S  as  the  portion  of  z  =  0  bounded 
by  4>,  we  have 

current  =  ff  =  _/[-(2=0)]  =  ./[Kj0ci<*  (23b) 

S  4>  <J> 

in  terms  of  the  elementary  vector  area  dS  of  S  and  the  elementary  vector  length 
dl  of  <t>.  The  use  of  X,  does  not  mean  that  this  current  is  flowing  in  the  z  =  0 
plane;  it  is  the  z-directed  current  flowing  into  or  out  of  that  plane  (in  the 
positive  z-direction) . 

These  definitions  of  voltage  and  current  can  readily  be  generalised  to  other 

two-conductor  systems.  For  one-conductor  systems,  however,  there  are  no 

"natural"  definitions  available,  and  it  is  necessary  to  appeal  to  arbitrary 

definitions  based  on  reflection  coefficients;  the  fundamental  relation  here  is 

that  giving  the  reflection  coefficient  T  of  an  admittance  Y  normalised  to  a 

characteristic  admittance  Y  (see  [  1]  ,  [  3) ) 

o 

Y_  _  1  -  r 

y  l  +  r 

o 

This  relation  can  be  used  in  conjunction  with  eqns  (21)  and  (22)  to  obtain  Y, 

by  the  method  described  in  l  1]  ;  but  this  method  involves  several  complications, 

in  particular  the  introduction  of  the  concepts  of  reaction  and  complex  power, 
and  the  implied  assumption  that  the  mode  constants  a.  ,  6.  are  all  real  (or  at 
any  rate  can  be  made  so  by  multiplying  them  all  by  the  same  complex  number). 

We  will  avoid  these  problems  by  continuing  the  development  independently  of 
[  1]  (an  example  of  the  connection  between  the  two  treatments  will  be  given  later). 

9  THE  INTEGRAL  RELATIONS 

At  this  point  we  introduce  four  essential  integral  relations  between  the 
mode  functions.  It  is  well  known,  and  is  proved  in  [  1]  ,  [2]  and  [3],  that  the 


mode  functions  appearing  in  such  equations  as  (17)  are  orthogonal  with  respect 
to  the  cross-sectional  area  of  the  generalised  cylindrical  system  to  which  they 
relate;  that  is,  if  fj(p,$)  and  f ^ ( P » <{>)  are  two  such  transverse  mode  functions, 
then 


transverse 


If 


f1(p.4>)f2(p*4,)dA 


cross-section 


0  unless  f^  = 


In  the  present  case,  where  dA  is  pdpdip  and  the  functions  are  all  independent 
of  this  suggests  •'hat  integrals  of  the  form 


f 


f1(p)f2(p)pdp 


will  be  important.  Regarding  the  TEM  mode  function  as  being  in  a  class  of  its 
own,  we  obtain  four  such  fundamental  integrals,  all  of  which  can  be  put  in 
closed  form  using  (  4]  (similar  results  also  exist  in  other  problems  with 
generalised  cylinders,  but  they  cannot  usually  be  expressed  so  conveniently; 
see  [  21 ) .  The  relations  are 


(24a) 


f 


Z“(0p)pdp 


!  P2(z2(eP)  -  zo(eP)z2(0P)) 


(24b) 


(ep)pdp 


-  i  zo(9o) 


(24c) 


(6p)Z^($p)pdp 

(e*») 


(4>pz.(eP)z  Up)  -  dpz,(4>p)zn(0p)) 

10  to 


(24d) 


Eqns  (24c)  and  (24d)  are  the  orthogonality  relations  properly  so  called,  and 
when  suitable  limits  are  inserted  and  eqn  (14)  used  it  can  be  seen  that  the 
mode  functions  are  indeed  orthogonal.  (The  definitions  of  the  Z  functions  are 
taken  from  eqns  (13).)  Eqn  (24d)  is  valid  even  if  6  and  $  belong  to  different 
families  (one  a  k , .  and  the  other  a  k^)  although  orthogonality  does  not  then 
hold  good. 
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10  CURRENT  AND  VOLTAGE  AT  THE  PLANE  OF  DISCONTINUITY 


From  eqns  (14),  (21),  (22),  (23)  and  (24c)  we  can  now  determine  the 
voltages  across  the  discontinuity  plane,  and  the  currents  flowing  into  it  and 
out  of  it,  on  the  A  and  B  sides.  In  an  obvious  notation  we  have 


R  R 

f  lX]dp  =  /  p  ao(1  +  rA)dp  =  a0(1  +  rA)ln(R/rA) 


(25a) 


and 


R  R 

■  /  [*»]*>’/ +  rB)dl> 


6o(l  +  rB)ln(R/rfl) 


(25b) 


Also 


R  R  j 

r  r 

00 

II 

< 

> 

1  7  “o(1  *  rA)d0  ‘  /  , 

—  a  (1  +  r  )  +  N  '  a.Z  (k  .p) 
p  o  A  /  >  l  1  Ai 

r 

A  rA  ' 

i=l 

dp 


(using  eqns  (24c)  and  (14),  and  inverting  the  order  of  integration  and  summation 


in  Fourier  anal 

ysis) 

R 

R 

Ri 

00 

■  /  ‘ 

I  ^  " 

/ 1 

7  60(1  8jzt(kBj0> 

rA 

rB 

j-1 

dp 


from  eqn  (21) 


K 

/ 


7  s  (i  +  r  )dP 

P  O  D 


from  eqns  (24c)  and  (14) 


=  Vg  =  V  (say) 


(25c) 
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Again 


j  K]  pd*  ■  j  ^ », 
/'  W-’/i 


—  «0<i  -  rA)pd<. 


—  o  (i  -  r.) 

n.  o  A 

A 


(26a) 


A  »A-  1  +  r  )pd4> 

p  o  o 


-  12.  6  (1  -  r  ) 

nB  °  B 


(26b) 

Since  V  =  Vfl,  it  is  clear  that  the  equivalent  circuit  representing  the 
discontinuity  must  consist  simply  of  a  shunt  admittance;  any  other  components 
would  introduce  a  voltage  di  ontinuity  across  z  ■  0.  (In  the  truncated  system 
(Fig  la)  Vg  has  no  meaning;  b  it  in  this  case  there  is  no  free  propagation 
beyond  2=0  and  the  equivalent  circuit  of  the  whole  system  is  terminated  there, 
so  we  again  arrive  at  a  shunt  admittance  as  the  representation  of  the  dis¬ 
continuity.)  A  current  IA  arrives  at  z  =  0  from  the  A  side  and  a  current  1^ 
leaves  2=0  from  the  B  side  [since  I.  and  I_  are  both  flowing  in  the  positive 
z-direction,  by  eqn  (23b)).  ihe  difference  of  these  two  currents  must  be 
considered  to  be  driven  through  the  shunt  admittance,  Y,  by  the  voltage 


VB  *  V;  so 
YV  =  1  - 


(27) 


This  equation  also  applies  to  the  system  of  Fig  la,  since  then  *  0 

(as  remarked  after  eqn  (17a))  and  so  from  eqn  (26b)  I  =  0;  eqn  (27)  then 

o 

st  ites  that  the  entire  curr  ;nt  IA  arriving  at  z  =  0  from  the  A-side  is  diverted 
th  ough  the  shunt  admittanc  Y,  which  is  clearly  true  since  in  this  system 
th  i  equivalent  circuit  is  a  terminating  impedance  at  z  =  0.  In  this  case  the 
sinplified  form  oT  eqn  (27)  YV  =  IA>  may  be  reduced  using  eqns  (25)  and  (26) 
to 


which  illustrates  the  ;ener  l  relation  from  [  ll  stated  previously, 

Y_  _  1  -  r 

Y  1  +  r 
o 


with  the  characteristu  admittance  Yq  assuming  is  transmission-line  value. 
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If  we  substitute  into  eqn  (27)  from  eqns  (25)  and  (26),  we  obtain 


(28) 


using  eqn  (22)  and  introducing  the  dummy  variable  p'  to  avoid  confusion  with 
p  (which  may  assume  any  value  in  the  range  r  <  p  <  R;  for  p  values  within 

A 

this  range  the  right-hand  side  of  eqn  (28)  is  actually  a  constant,  although 
it  looks  like  a  function  of  p).  The  integrations  in  this  equation  and  in  the 
equations  which  follow  extend  over  the  gap  between  the  conductors  at  z  =  0. 

In  the  present  case  the  dimensions  oi  this  gap  coincide  with  those  of  one  of 
the  pairs  of  cylinders,  but  this  need  not  be  the  case  (for  instance,  there 
might  be  a  diaphragm  at  z  =  0;  or  the  cylinders  might  be  stepped  so  that  their 
radii  both  increased,  which  would  create  a  gap  smaller  than  the  spacing  between 
the  cylind  ;rs  on  either  s  de,  as  in  Fig  le) .  The  integrations  extend  formally 
over  the  area  of  the  gap,  and  to  apply  the  method  of  this  memorandum  in  general 
it  is  necessary  to  replace  the  single  integrals  in  eqns  (25) ,  (26)  and  (28)  by 
double  integrals;  the  formulation  in  this  case  is  discussed  in  [1].  However, 
the  rotational  symmetry  of  the  systems  represented  in  Figs  1  allows  us  to  deal 
with  them  using  the  simpler  single  integrals. 

11  EXPRESSION  FOR  THE  SHUNT  ADJUTTANCE 


We  now  express  the  quantities  ou  and  0.  in  terms  of 
we  have 


From  eqn  (21) 


/R 

42l(k#i»)p«  •  /  i.0(l  ‘  V  *  Z  VZl(kAC 
r.  r.  i'=l 


p) 


Zl^kAip^pdp 


(introducing  the  dummy  suffix  i'  to  avoid  confusion  with  i) 

•  /  “ o(1  *  rA,2l(kAi»,d»  *  £  /  VVaC^V11 

r.  i'-l  V. 


p)p  dp 
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(interchanging  the  order  of  integration  and  summation) 


•.«  *  V  ['  W’]  * 


Ev 


i'=l 

i'*i 


(kAipZl(kAi'p)Zo(kAiP>  ‘  kAi*pZl(kAiP)Zo(kAi'P)) 


~i  R 


<kAi-  -  kAl> 


R 

+a i  ^  ZZ(kAip)pdp 


(using  eqns  (24c)  and  (24d)) 


«i/ «? 


(kAip)PdP 


(29a) 


(using  eqn  (14)) 


This  is  an  example  of  the  orthogonality  property  of  the  mode  functions.  The 
integral  on  the  last  line  of  eqn  (29a)  can  also  be  written  in  closed  form 
using  eqn  (24b),  but  we  shall  leave  it  as  it  is  for  the  time  being.  We  can 
also  write,  using  eqn  (21), 


£pZl(kBjp)pdp 


R 

f  £pZl(kBjp)pdp 
rB 


(29b) 


In  terms  of  our  earlier  discussion  of  these  integrations,  the  last  equation 
transforms  an  implicit  double  integral  over  the  gap  at  z  *  0  to  one  over  the 
whole  cross-section  there.  Since  the  cross-section  consists  entirely  of  the 
gap  and  a  flat  section  (or  sections)  of  perfect  conductor,  such  a  transformation 
is  always  possible,  even  when  the  system  is  not  rotationally-symmetrical  (see 
[  1] ) .  In  systems  like  that  of  Fig  le,  it  is  necessary  to  transform  in  this  way 
on  both  sides  of  z  *  0,  and  we  then  have  to  introduce  a  third  set  of  orthogonal 
functions  (orthogonal  over  the  gap,  instead  of  over  one  of  the  generalised 


cylindrical  cross-sections).  In  the  present  problem,  fortunately,  this 
complication  is  unnecessary. 

From  eqns  (29b)  and  (21)  ,  we  now  obtain 

R  R  | 

/  Spzi(kBjp)|Id0 '  /  h  6o(1  *  V  *  S  VWP)  vs/*04’ 


>j  J 


Bj  /  Zjftypjpdp 


on  repeating  the  steps  in  the  derivation  of  eqn  (29a).  Substituting  eqns  (29) 
into  eqn  (28)  gives,  using  the  dummy  p'  as  before. 


K 

/f  ' Zl^kAiP  *  ^  ^  * d(3 


y  zj(kAip')p’dP’ 


Z!(kAip)  + 


R 

-  J*’'1' 

mtj: 

J  zl<k 


(kgj  p  ' )  p '  dp ' 


Zj(kBjp)  •  (30) 


(kgj  p  ' )  p '  dp ' 


Finally,  if  we  multiply  this  equation  throughout  by  £  p  and  integrate  across 
the  gap  (which  in  our  reduced  terms  is  equivalent  to  multiplying  by  £  and 
integrating  over  the  area  of  the  gap)  we  obtain,  on  interchanging  the  order 
of  integration  and  summation, 


'  R  \  2  ® 

/'■*■  -mt! 


(f  <?pZl(kAiD)pdD) 

j-A^A  / 


/  Z‘<kA‘ 


P)  pdp 
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12  STATION ARITY  AND  OTHER  VARIATIONAL  PROPERTIES 


(31) 


Eqn  (31)  is  a  special  case  of  a  general  formula  derived  in  [  1]  ,  a  so-called 
variational  formula.  The  properties  of  variational  formulae  are  discussed  in 
[  1]  and  { 2] ,  where  it  is  shown  that  they  have  the  remarkable  property  of 
stationari ty.  In  our  case  this  means  that,  if  the  true  field  ^  which  appears 
in  eqn  (30)  is  replaced  by  an  arbitrary  function  £  ,  then,  provided  the 
arbitrary  function  (like  the  true  field)  satisfies  the  proper  boundary  conditions 
at  the  edges  of  the  gap. 


— =  0  when  £  =  £ 

»*.  c  0 

A  proof  of  this  property  is  most  conveniently  given  using  the  development 
of  the  calculus  of  variations  by  Courant  and  Hilbert  [5].  They  show  that  a 
double  integral  of  the  form 


b  b 


K(p.p')  <?Q(p)  £p(p')dpdp' 


a  a 


(32a) 


where  a  and  b  are  constants,  and  K(p,p')  is  a  fixed  function  satisfying  con¬ 
ventional  conditions  of  continuity  and  differentiability  and  possessing  the 
symmetry  property  K(p,p')  ■  K(p',p),  will  be  stationary  with  respect  to  changes 
in  the  functional  form  of  if  its  functional  gradient,  which  is 

b 

2  J  K(p ,p ' )  (p ' )dp '  ,  (32b) 

a 


vanishes  for  all  values  of  p  in  the  range  a  <  p  <  b.  The  squared  integrals  in 
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eqn  (31)  can  be  written  as  double  integrals  of  the  form  (32a);  freely  inter¬ 
changing  the  orders  of . integration  and  suinnation,  and  replacing  p  by  p'  where 
appropriate,  we  find  thai.  eqn  (31)  and  the  equation  obtained  by  setting 
expression  (32a)  equal  to  zero  can  be  identified  with  each  other  if  we  choose 


a  “  r .  ,  b  -  R  and 
A 


K(p.p’) 


Zl(kAip)pdp 


j  V^VaI 


P')PP'  + 


Zl^kBjP^pdp)  Z1(kBjp)Z1(kBjP,)pp'  “ 


Y_ 

2v 


(32c) 


which  has  the  required  symmetry  property  under  interchange  of  p  and  p'.  It 
them  turns  out  that  the  fi nctional  gradient,  expression  (32b),  vanishes  when 
£  ~  £p>  because  of  eqn  (30). 

Now,  suppose  we  try  to  find  an  approximate  value  of  Y  b\  using  an  approxi¬ 
mate  function  £  to  represent  £  in  eqn  (31).  If  the  error  in  the  approximation 

~  ^  P 

£  is  of  the  first  order  of  small  quantities,  the  error  in  the  computed 
approximate  Y  will  only  be  of  the  second  order  of  smallness,  because  of  the 

/K  A 

vanishing  partial  derivative  with  respect  to  £^  at  £^  «  .  This  is  an 

important  consequence  of  stationarity;  but  still  more  im>ortant  is  the  fact 
that  stationarity  gives  us  a  simple  dm  t hod  of  constructing  good  approximations 
£  •  if  we  imagine  generalised  cylinders  to  be  obtained  by  translating  along 
the  z-axis  the  boundary  curves  of  the  gap  in  the  plane  z  =  0,  then  the  mode 
functions  appropriate  to  these  cylinders  satisfy  the  proper  boundary  conditions 
in  the  plane  of  the  gap,  and  so  any  linear  combination  of  these  mode  functions 

is  a  possible  approximation  £  .  If  we  choose  a  finite  linear  combination,  use 

*  P 

it  as  an  approximation^  to  represent  £  in  eqn  (31),  and  adjust  all  the 
linear  multipliers  tc  satisfy  the  stationarity  condition  (3Y/>£p)  *  0,  we 
shall  obtain  a  second-order  approximation  to  the  value  of  Y.  By  using  more  and 
more  mode  functions  in  linear  combination,  we  can  obtain  a  sequence  of  approxi¬ 
mations  to  the  value  of  Y.  Since  the  mode  functions  form  a  niathematically 


"complete"  set  (see  l  2]  or  (  5] ) ,  we  can  represent  the  true  as  accurately 
as  we  please  (at  least  in  the  mean,  which  is  all  that  matters)  by  taking  a 
su£cicient  number  of  them  in  a  suitable  linear  combination;  so,  if  the 
sequence  of  approximations  tends  to  a  limit  (which  we  shall  show  to  be  the 
case),  the  limit  will  be  the  true  value  of  Y.  This  is  the  well-known  Ritz 
procedure  (discussed  in  (  11  ,  [2]  and  [  5] ) . 

13  APPROXIMATIONS  BY  THE  RITZ  PROCEDURE 

We  now  apply  the  procedure  just  described.  Let  £  in  eqn  (31)  be  replaced 
~  P 

by  an  approximation  £  of  the  form 


£  =  —  a  +  /  ^  k  .  ,a.  ,Z,  (k  . ,( 

p  P  o  Ai  l  1  Ai 


since  the  generalised  cylinders  constructed  by  translating  along  the  z-axis 
the  gap  in  the  plane  z  =  0  are,  in  this  case,  just  the  same  as  the  cylinders 
on  the  A-side  (if  they  were  different  from  the  cylinders  on  the  A  and  B  sides, 
as  in  the  system  of  Fig  le,  it  would  be  necessary  to  introduce  a  third  set  of 
orthogonal  mode  functions,  but  in  the  present  problem  this  complication  does 
not  arise).  The  factor  ,  in  each  term  of  the  sum  in  eqn  (33)  is  inserted 

for  future  convenience,  and  N  is  an  unspecified  positive  integer.  We  then  have, 
from  eqns  (24)  and  (14) , 


K 

/  zi;k«0) 


I  -  rIZl(kAirA» 


(34a) 


/ 


(34b) 


a  In  (R/r  ) 
o  A 


(34c) 
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R 

/  4Zl(kAiP)pdp 


T  kAi;i(R2z?(kAiR) 


riZl(kAirA)) 


(i  <  N) 


(34d) 


=  0 


(i  >  N) 


f  ~  -  Z  (k  .r  )  JLv 

J  €< pZl(kBjp-pdp  =  -  +  2 


Bj 


2  ,2  J  kBjrAZl(kAi'rA)Zo(kBjrA) 


i'=l  \kBj  ~  kAi' 


(34e) 


On  making  use  of  eqns  (33)  and  (34)  in  eqn  (31),  we  obtain  an  approximation 
to  Y,  Y,  given  by 


_Y _ 

2  it  j  oi 


In  (R/rA)j2 


t  ft) 


(kA.R) 


ri?XirA> 


■) 


.  2  '2 
k.  . a. 
At  i 


which  becomes,  on  interchanging  the  order  of  summation  and  making  some  convenient 
changes  in  the  dummy  variables, 


(a~  in  (R/rA))2  - 


2irj  oj  '  O 


j  =  l 


4Z  ik  .r  ) 
o  Bj  A 


(k0jRZi  (kgjH))  ~  (kBjrBZl^kBjrB^ 


Z]  “i  [iTT  (<kAiRZl<kAi 
i-1 


AiR»2  -  <kAikAZl(kAirA»^l  * 


)] 
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+  2a 


o  /  j  1 


i=l 


(kAirAZl(kAirA>) 


kii 


~  k?,  /TV 


j=l  \  Bj  Ai 


*Zo<kBjrA> 


k(kBjR2i(kBjR))  ~  (kBjrBZl^kBjrB')') 


N  N 


Z  Z  Vi- 

i=l  i’=l 


^kAirAZl^kAirA^  ^kAi,rAZl^kAi,rA^  ^  iZ 


Ji 


2  2 

.  ,  \kn.  -  k". 
j=l  \  Bj  kx. 


'B 


^  Bj/  l  ^Bj^l^Bj^  ”  ^kBjrBZl^kBjrB^ 


(35) 


14  SIMPLIFICATION  OF  THE  APPROXIMATE  EXPRESSION 

Eqn  (35)  admits  a  number  of  simplifications.  In  the  first  place,  it  is 
clear  that  only  the  ratios  of  the  a’s  are  significant,  since  the  equation  is 
homogeneous  in  them.  It  is  convenient  to  write 


“i^AiVl^Ai^  =  "  Vo 


(36a) 


where  the  a^  are  the  new  unknown  multipliers. 

Secondly,  if  £  denotes  one  of  the  quantities  (k  .r  ),  (kD.r0),  (k..R),  (k„.R), 

Al  A  DJ  o  Al  D J 

we  have  from  eons  (12)  and  (13) 

J  (O 

ZjU)  =  J^s)  -  yTo-  yiU) 
o 

But  for  any  complex  number  w,  it  is  true  that 

J, (w) Y  (w)  -  J  (w)Y  (w)  *  J  (w)Y’(w)  -  Y  (w)J’(w)  »  — 

io  ol  oo  oo  7TW 

(after  (2]  or  [  4]  ) 

provided  w  ¥=  0,  where  the  prime  denotes  differentiation;  so 


26 


(36b) 


ZX(D  =  2/(ttCYq(0) 

(£  *  0,  because  the  k^,  the  kgj  and  rA>  rfi  and  R  are  all  strictly  positive). 

.  .  ‘2 

We  therefore  have,  after  dividing  by  aQ, 


2irjw 


N 


-‘£-i 


i=l 

N  N 


EE 

i=l  i'=l 


(from  eqn  (36a)) 


E 

j=i 


•^foniV 


2TBj/(VkBjR>r2  '  (Yo<kBjV> 


N 

E  4 

i=i 


eA  //VkAirA) 


^AiWY^k^R) 


N 


-  2 


E\ 

i  =  ] 

N  N 

EE 

i=l  i ' =1 


(from  eqn  (36b) ) 
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Again,  from  eqns  (12)  and  (13)  we  have 


Zo(kBjrA)  ”  Jo(kBjrA)  "  Y  (k  "1 R)  Yo*kBjrA) 

O  BJ 


so  if  we  define 


x. 

J 


(J  (k  . r  )Y  (k  . R)  -  J  (k„ . R) Y  (k_.rA))‘ 
o  Bj  A  o  Bj  o  Bj  o  Bj  A 


rBj 


/VSj!L\ 

VWb  V 


(38a) 


then 


x.  * 
J 


V<Yo(V>r2  •  <’fo'kBi,:B))'2) 


(38b) 


If  we  also  introduce  the  relative  permeabilities  and  permittivities  p  , 
__  _ _  A 

Ug ,  eA»  Eg,  the  absolute  permeability  and  permittivity  of  free  space  p^,  eQ, 
the  velocity  of  light  in  free  space  c  and  the  well-known  relation  ([  1]  ,  [2]) 


pec 
o  o 


we  obtain,  using  eqn  (6a), 


2 - .  2 

“  PA  CA  7  C 


(38c) 


with  a  similar  relation  for  k„,  and 

D 


eB  =  eB  7  <PoC  ) 


(38d) 


Introducing  eqns  (38b)  and  (38d)  into  eqn  (37),  along  with  the  further 
definitions 


j-1 


(38e) 
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gives  us 


(38f) 


(38g) 


(38h) 


Y_ 

jw 


u  c2(ln  (R/r  ))2 
o  A 


E2 

a.  s. 
i  i 


i=  1 


N  N  N 

a.  t .  + 
l  l 

i=l  i=l  i ' =1 


(39) 


15  STATIONARY  VALUE  OF  THE  APPROXIMATE  ADMITTANCE 


Applying  the  Ritz  procedure  to  determine  the  stationary  value  of  Y  in 
eqn  (39)  is  now  just  a  matter  of  taking  its  partial  derivatives  with  respect  to 
each  of  the  in  turn  and  setting  all  the  resulting  expressions  equal  to  zero. 
Taking  note  of  the  fact  that  u.  .,  =  u.,  .  (eqn  (38h))  it  is  easy  to  see  that 


3 

3a* 


■1 


N 

*  £ 

i=i 


£ 

L-l 


a.  t. 
l  i 


N  N 

££ 

i=l  i'=l 


a .  a.  .  u . 


i.i 


-  2t, 


+  2Vm 


+  2 


i=l 


and  consequently  the  values  of  the  a^  which  make  Y  stationary  are  given  by 
solving  the  simultaneous  linear  equations 


(s 


+  Ui,i) 


e,i 


a . 
i 


a  =  i,  2 . n) 


(40) 


i-1 
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where  a  superscript  N  has  been  introduced  on  the  a^  to  stress  that  their  best 
values  depend  on  N.  Substituting  eqns  (40)  into  eqn  (39),  and  writing 
for  the  stationary  value  of  Y  (which  depends  on  N,  because  the  best  values  of 
the  a.  do  so)  we  obtain 


V  c2(ln  (R/r.))2 
O  A 


(41) 


because  multiplying  each  of  eqns  (40)  by  the  corresponding  a^  and  then  adding 
all  the  results  together  gives 


E' 


2 

V* 


4=1 


N  N 

♦EE'-. 

4=1  i=l 


a .  u.  . 
l  4,i 


a.t. 
4  4 


4=1 


It  can  be  seen  that  the  solutions  of  eqns  (40)  are  real  numbers  (not 

complex  ones),  because  the  definitions  in  eqns  (38)  show  that  the  known 

quantities  in  eqns  (40)  are  all  real.  It  is  then  clear,  from  eqns  (36a)  and 

(33),  that  the  approximate  field  £  is  real,  apart  possibly  from  the  connon 

factor  a  (this  implies  that  the  fields  due  to  the  individual  modes  are  all  in 
o 

phase  with  each  other,  which  is  a  quite  reasonable  physical  condition  for 
stationarity) .  Since  £  can  approach  &  as  closely  as  we  please,  and  since  the 

a  P  P  A. 

factor  aQ  is  of  no  significance  in  any  expression  which  is  homogeneous  in  ^ , 

it  then  follows  that  6^  in  eqn  (31)  can  be  taken  as  a  pure  real,  and  consequently 

that  Y/(jw)  is  a  pure  real.  Moreover,  this  pure  real  is  positive,  because  when 

e<-n  (31)  is  rewritten  as  an  expression  for  Y/(jw)  it  becomes  apparent  that 

every  quantity  in  the  expression  is  a  perfect  square  or  otherwise  known  to  be 

positive  (including  the  integrals  in  the  denominators,  which  have  non-negative 

integrands).  It  is  clearly  also  true,  by  the  same  argument,  that  every  YN/(jw) 

is  a  positive  real.  It  follows,  therefore,  that  the  quadratic  form  in  the  a^ 

for  Y,  eqn  (39),  must  be  positive-definite,  and  that  it  attains  an  absolute 

minimum  when  the  best  values  given  by  eqn  (40)  are  inserted  in  it.  If  N  is 

increased  to  (N  +  1) ,  then  in  terms  of  the  new  approximate  &  the  old  stationary 

P  N 

field  is  no  longer  stationary,  and  when  it  is  made  so  by  replacing  the  a^  by 
the  ^N+^a.  it  attains  a  new  absolute  minimum  which  must  be  lower  than  the  old 

l 

one.  Accordingly  the  Y^/(ju)  form  a  monotonic  decreasing  real  sequence  in  N, 
which  is  bounded  below  (since  every  member  of  it  is  positive);  and  so,  since 
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£  can  be  made  to  approach  as  closely  as  we  please  by  taking  N  large  enough, 
it  follows  that,  as  N  tends  to  infinity,  YN/(jw)  tends  from  above  to  the  value 
given  by  eqn  (31)  when  ,  whicli  is  Y/(joi). 

16  THE  DISCONTINUITY  CAPACITANCE 

We  can  now  complete  the  formal  processes  by  providing  an  interpretation  of 

Y.  We  have  already  shown  that  Y/(jw)  is  a  positive  real,  and  since  the 

admittance  of  a  capacitance  is  (jwC)  when  the  time  variation  of  the  fields  is 

as  exp(jwt)  (as  assumed  at  the  beginning)  it  is  convenient  to  interpret  Y  as  the 

admittance  of  a  capacitance  Y/(jw).  This  becomes  not  only  convenient  but  also 

natural  when  we  note,  from  eqns  (10),  (15c),  (38),  (39)  and  (40),  that  Y/(ju>) 

2  2 

depends  only  slightly  on  the  frequency  (through  the  influence  of  kA  and  kg 

on  the  yA£  and  Ygj)  when  the  frequency  is  well  below  the  upper  critical 

frequency  defined  earlier,  and  that  Y / ( j co)  attains  a  definite  (positive)  value 

when  the  frequency  is  zero  (from  eqns  (15c)  and  (38c)  we  then  have  simply 

y..  =  k..,  y„.  =  k„ .  and  there  are  no  numerical  difficulties).  Also,  multiplying 
Ai  Ai  Bj  Bj 

e  and  e  simultaneously  by  a  given  factor  causes  the  zero-frequency  value  of 
A  B 

Y/(jw)  to  increase  by  the  same  given  factor,  a  result  which  remains  approximately 
true  at  non-zero  frequencies  well  below  the  upper  critical  frequency.  These 
are  all  characteristics  of  a  slightly  imperfect  capacitance;  so  we  conclude  that 
the  equivalent  circuit  of  a  system  like  that  of  Fig  la  or  Fig  lb  should  be  a 
shunt  capacitance,  having  a  value  given  by  (Y^/ijui))  with  Y^  defined  by 

eqn  (41).  A  similar  conclusion  holds  good  for  the  other  systems  of  Figs  1. 

Most  of  the  theory  of  the  last  few  sections  can  be  generalised  to  more 
complicated  systems,  as  discussed  in  (  1]  ;  in  particular  it  is  nearly  always 
possible  to  establish  a  variational  equation  like  eqn  (31),  although  the 
generalised  variational  equation  will  not  necessarily  lead  to  a  positive- 
definite  quadratic  form.  If  no  such  form  can  be  found,  the  smoothness  and 
speed  of  convergence  to  the  true  value  will  be  impaired,  but  convergence  will 
still  occur. 


The  special  case  represented  by  Fig  la  can  be  treated  by  allowing  r^  to 
tend  to  0  in  the  formulae  derived  above;  eqn  (10b)  reduces  to  eqn  (20a),  and 

the  functions  z0/q/2^kBjP^  reduce  t0  Jo/l/2^kBjP^  ’  froro  e<ins  and  and 

the  discussion  preceding  them  (note  that  Y^k^rg)  tends  to  infinity  as  rg 
tends  to  0).  Eqn  (38a)  then  assumes  the  special  form 
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<Jo<kBjrA)To<kBt*»' 


No  other  changes  are  necessary. 

17  PRACTICAL  DETAILS  OF  THE  CALCULATION 

It  only  remains  to  say  how  eqns  (10),  (15c),  (20a),  (38),  (40),  (41)  and 
(42)  should  be  used  in  practice  to  compute  the  shunt  discontinuity  capacitance. 

The  infinite  sums  in  eqns  (38)  are  naturally  truncated  to  large  finite 

upper  limits.  If  the  last  term  retained  is  the  Mt\  the  error  committed  by 

-2 

truncation  is  of  the  order  of  M  ,  as  can  be  shown  by  the  following  argument. 
Computing  some  of  the  sums  for  typical  cases  shows  that  the  most  important  term 
in  the  brace  brackets  in  eqn  (41)  is  the  leading  term  q,  as  we  expect  from  the 
fact  that  when  N  =  0  the  expression  in  the  brace  brackets  actually  reduces  to 
q.  It  follows  from  eqn  (38e)  that  we  can  estimate  the  truncation  error  by 
simply  considering  eqns  (38a)  and  (42) .  The  discussion  in  [  2]  shows  that  the 
large  values  of  kfij  are  of  the  asymptotic  form  j7r/(R  -  rB>  (this  formula 
actually  holds  fairly  well  also  for  small  j,  and  was  used  to  derive  the  approxi¬ 
mate  formulae  for  the  critical  frequencies).  From  [  2J  or  [  4] ,  J  (z)  and  Y  (z) 

.  o  o 

fall  off  like  z  i  for  large  z,  and  from  eqn  (15c)  yn .  is  asymptotic  to  k  .  for 

DJ  Dj 

large  j.  It  follows  from  eqns  (38a)  and  (42)  that  the  x.  fall  off  like 

-3  .  J  , 

j  for  large  j,  and  the  relative  error  committed  by  stopping  at  some  large  j 

equal  to  M  is  therefore  of  the  order  of 


Ei 

j=M+l 


or  of  the  order  of  M  ,  using  integrals  to  approximate  the  sums.  The  error  can 
therefore  be  reduced  to  something  of  the  order  of  10  ^  by  taking  M  ■  1000,  as 
was  done  by  Somlo  in  the  program  referred  to  in  his  paper  [  6] ;  this  is  acceptably 
small,  both  in  its  own  right  and  in  comparison  with  Somlo's  estimate  of  the 
accuracy  of  the  whole  calculation  (which  appears  to  depend  mainly  on  the  extra¬ 
polation  to  infinite  N) .  We  therefore  cut  off  all  the  infinite  sums  after  the 
first  1000  terms. 
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The  Bessel  functions  can  be  evaluated  by  many  standard  methods  and  we  shall 
not  discuss  them  in  detail.  It  is  worth  noting,  however,  that  since  the 
contribution  of  the  j1"*1  term  to  each  sum  falls  off  in  importance  like  j  ,  it  is 
acceptable  to  use  methods  which  give  progressively  worse  accuracy  in  Jq(z)  and 

Y  (z)  as  z  increases  (such  as  the  fast  and  fairly  simple  method  based  on  a 

0  .  -1 
Chebyshev  series  in  z  for  small  z  and  an  asymptotic  expansion  in  z  for  large 

z,  butted  together) . 


The  quantities  k..  and  kn .  can  be  computed  from  eqns  (10)  and  (20a),  using 
Al  Dj 

standard  methods  of  finding  zeros.  Since  it  is  necessary  to  find  every  zero 
in  a  given  range  once  and  only  once,  it  is  desirable  to  use  a  method  which  takes 
an  upper  and  a  lower  limit  and  finds  a  zero  between  them,  using  Bessel  function 
theory  to  make  sure  that  there  is  one  and  only  one  zero  between  the  given 
limits.  Suitable  limits  may  be  obtained  from  the  asymptotic  formula  quoted 
above,  and  its  degenerate  form  for  the  system  of  Fig  la  which  is  obtainable  to 
sufficient  accuracy  by  putting  r  =  0  (see  [  2]  and  [  4] ) .  We  then  find  that  the 
solutions  of  eqns  (10b)  and  (20a)  each  satisfy  the  relation 


(R^)  -  °-3)  <  kBj 


< 


(R  -  rB) 


(j  +  0.3) 


(43) 


where  the  (somewhat  arbitrary)  choice  of  ±0.3  is  wide  enough  to  include  one 
zero  and  narrow  enough  to  include  only  one  (and  that  one  the  jth>.  (From  [4]  , 
there  are  no  multiple  zeros.)  In  the  same  way  we  can  set  limits  to  k^  for  use 
with  eqn  (10a). 

The  solution  of  eqns  (40)  for  a  given  N  can  also  be  obtained  by  standard 
methods.  Because  of  the  positive-definiteness  of  eqn  (39),  eqns  (40)  have  a 
positive-definite  coefficient  matrix.  This  is  easily  proved  by  writing 
a^  =  a^/a^.+^,  when  the  positive-definite  expression  in  the  brace  brackets  in 
eqn  (39)  assumes  the  form 
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The  homogeneous  quadratic  form  in  brackets  here  is  necessarily  also  positive- 
definite,  and  its  associated  matrix  B,  of  order  (N  +  1)  x  (N  +  1) ,  is  therefore 
positive-definite  too  and  has  elements  b.  . ,  satisfying  the  conditions 

^  i  ^ 

bi  i’  =  ^ ,  i  an<*  (fro®  expression  (44a)) 


N+l  N+l 


ZEai  bi,i-  v 


Vi  * + 


Vr2..  .2yrt.r;  + 

JLmJ  1  1  /  1  1  1  N+l 


i*l  i'“l 


N  N 


IL‘i  \.v  V 


(44b) 


i=l  i’=l 


Equating  corresponding  terms  and  taking  account  of  the  fact  that  u^  =  u., 
we  have 


From  Sylvester's  theorem  (see  Faddeev  and  Faddeeva  [  7] )  all  the  leading 
sub-matrices  of  a  positive-definite  matrix  are  also  positive-definite,  in 
particular  the  N**1  such  sub-matrix  of  B,  which  consists  of  the  elements  b^ 

(i  <  N,  i'  <  N) ;  and  from  eqns  (44c),  this  sub-matrix  is  precisely  the 
coefficient  matrix  of  eqns  (40).  Special  algorithms  exist  for  solving  a  system 
of  equations  with  a  positive-definite  coefficient  matrix,  which  reduce  by  half 
the  computational  effort  required;  however,  none  of  these  algorithms  were  used 
in  programming  the  theory  developed  above,  for  reasons  which  are  explained 
later.  The  positive-definiteness  has  another  and  much  more  important  con¬ 
sequence:  it  guarantees  that  eqns  (40)  always  possess  a  unique  solution  (see 
[  7] ) .  This  solution  is  not  necessarily  we 11 -determined  (eqns  (38)  show  that  if 
any  one  of  the  k„ .  is  nearly  equal  to  one  of  the  k, . ,  some  of  the  u.  .,  and  t. 
will  be  very  large  and  not  accurately  calculable),  but  in  practice  i  is  always 
small  enough  for  this  problem  not  to  arise.  The  author's  program  includes 
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a  check  that  no  difference  (k.  .  -  k„.)  is  less  than  10  of  the  k  used 

Ai  oj  Ai 

in  computing  it,  and  this  check  is  always  satisfied  in  the  numerical  cases  so 
far  treated,  which  restricts  the  loss  of  significant-figure  accuracy  to  an 
acceptable  three  figures. 

With  the  information  in  this  section  it  is  a  straight-forward  matter  to 
construct  a  computer  program  which  will  calculate  Y^/(jco)  for  any  given  N, 
using  as  d?ta  the  values  of  w,  pA>  u  ,  eA,  e^,  R,  rA>  rfi,  and  c.  The 
extrapolation  to  infinite  N  will  now  be  discussed. 

18  EXTRAPOLATION  TO  INFINITY 

The  final  stage  of  the  computation  consists  of  determining  the  discontin¬ 
uity  capacitance  C,  Y/(jw),  from  the  values  of  Y^/ ( jou)  for  a  number  of  values 
of  N.  It  is  convenient  and  makes  good  mathematical  sense  to  choose 
consecutive  values  of  N;  in  [6],  for  instance,  the  values  28,  29,  ....,  40 
were  used.  The  theory  of  [6]  is  derived  from  Whinnery  et  al  [8];  the  final 
equations  arrived  at  by  these  authors  are  similar  to  eqns  (38)  to  (41),  but  no 
reasons  are  given  for  any  of  the  formal  manipulations  in  [  8] ,  the  variational 
property  is  not  mentioned,  and  the  results  are  rather  difficult  to  use  because 
they  have  not  been  simplified  as  fully  as  has  been  done  here.  (Some  of  the 
difficulties  of  use  are  mentioned  in  [  6] ,  and  these  difficulties  and  the 
inconveniently  concise  presentation  in  [8]  made  it  desirable  to  develop  the 
theory  systematically,  so  as  to  throw  light  on  the  numerical  behaviour  of 
the  approximate  solutions).  Because  the  program  corresponding  to  [6]  used  the 
theory  of  [  8] ,  which  is  similar  to  that  presented  above,  it  was  expected  that 
the  convergence  would  be  similar,  so  an  upper  limit  to  the  value  of  N  was  chosen 
at  40.  This  proved  to  be  satisfactory,  which  was  fortunate  because  the  time 
and  storage  demands  of  the  program  would  have  made  it  impossible  to  use  a 
substantially  higher  limit. 

In  (  6] ,  and  in  Bianco  et  al  (  9] ,  extrapolation  by  a  hyperbolic  fit 

'  -l 

(plotting  against  N  )  was  described,  but  the  program  corresponding  to  (6] 
(obtained  from  the  author  of  that  paper)  actually  used  an  implementation 
(no  215  of  the  ACM  Collected  Algorithms  [12))of  the  transformation  described 
by  Shanks  [  10]  t  and  this  more  powerful  method  was  employed  in  the  present  work. 
It  is  most  easily  applied  to  sequences  containing  (4n  +  1)  consecutive  elements, 
n  being  an  integer;  the  number  of  elements  in  the  set  28,  29 . .  40  is  of 
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this  form  with  n  *>  3,  which  suggests  that  Somlo  changed  his  procedure  while 
he  was  writing  his  paper.  (He  might  have  been  troubled  by  an  error  in  the 
printing  of  algorithm  215,  which  at  one  point  has  an  illegal  symbol  in  place 
of  a  semi-colon).  In  any  case,  we  are  provided  with  a  suitable  value  of  n. 

This  is  not  a  trivial  matter,  because  the  reduction  of  (4n  +  1)  elements  to  a 
single  final  value  actually  requires  n  applications  of  the  transformation, 
and  repeated  application  of  any  smoothing  algorithm  can  induce  instability  when 
the  differences  between  successive  smoothed  values  become  comparable  with 
the  errors  in  the  values  themselves;  examination  of  the  intermediate  values 
produced  in  the  author's  program  by  the  smoothing  algorithm  showed,  in  fact, 
that  taking  n  greater  than  3  would  have  caused  this  problem  to  show  itself. 

Values  of  C  calculated  using  N  »  28,  29,  . .  40  were  found  to  agree 

.  4 

with  the  tabulated  values  in  [6]  to  a  few  parts  in  10  .  The  author's  values 
were  consistently  higher,  which  suggests,  in  view  of  the  fact  that  the  true 
value  of  C  is  approached  from  above  as  N  increases,  that  the  values  in  [6] 
are  r ore  accurate.  However,  as  the  calculations  actually  performed  by  the  two 
programs  are  similar,  it  seems  that  the  estimated  accuracy  of  a  few  parts  in 
10^  for  the  values  in  [6]  may  not  be  attained.  For  practical  purposes  the 
discrepancy  is  too  small  to  matter.  It  is  also  worth  recording  that  agreement 
with  the  values  in  l  6]  to  a  few  parts  in  10  was  achieved  using  only  sets  of 

A 

Y^/  ( jui)  calculated  using  N  =  1,  2,  . .  5  with  one  application  of  the  e ^ 

transformation;  this  remarkably  high  accuracy  is  a  manifestation  of  the  fact 
that  the  are  second-order  approximations  to  the  true  value,  Y. 

19  OTHER  COMPUTATIONAL  POINTS 

It  can  be  seen  from  eqns  (40)  that  the  coefficient  matrix  in  the  linear 
equations  appropriate  to  any  given  N  contains,  as  a  leading  sub-matrix  of 
itself,  the  coefficient  matrix  appropriate  to  any  smaller  value  of  N.  One 
property  of  a  system  of  equations  having  a  positive-xlef inite  coefficient 
matrix  is  that  it  can  be  solved  accurately  (by  Gaussian  elimination  or  by 
resolving  the  coefficient  matrix  into  triangular  factors)  without  altering 
th*  order  in  which  the  equations  are  presented  (see  Wilkinson  [ 11]  ,  and  also 
[  7] )  so  that  solving  eqns  (40)  for  the  largest  relevant  N  does  most  of  the  work 
of  solving  them  for  all  smaller  values  of  N;  if  the  solution  is  undertaken  by 
triangular  factorisation,  for  instance,  the  triangular  factors  of  the  N'th 
leading  sub-matrix  of  the  matrix  for  the  largest  N  are  simply  the  N'th  leading 
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sub-matrices  of  the  triangular  factors  of  the  matrix  for  the  largest  N.  This 

fact  would  have  been  exploited  if  it  had  been  necessary  to  write  a  routine 

for  solving  simultaneous  linear  equations  instead  of  simply  calling  one  from 

the  computer's  library;  but  it  was  not  worthwhile  to  write  a  routine  specially 

because,  as  we  shall  now  see,  the  reduction  in  the  number  of  arithmetic 

operations  performed  during  the  whole  calculation  would  have  been  too  small 

to  be  useful.  From  [  7)  and  [ 11]  ,  an  ordinary  routine  for  solving  N  linear 

equations  in  N  unknowns  requires  about  N  /3  multiplications  and  the  same 

number  of  additions,  which  for  N  =  28,  29 . .  40  gives  about  2.10^ 

multiplications  in  total.  Examination  of  eqn  (38h)  shows  that  evaluating 

each  term  in  the  sum  giving  a  u.  requires  two  multiplications,  two 

» i 

divisions  and  two  subtractions,  in  addition  to  the  time  spent  in  determining 

2 

the  k.  . ,  k  .  and  x.;  there  are  N  / 2  distinct  u.  and  each  sum  contains 

AliJJ  J  1,1.. 

1000  terms;  so  for  N  =  40  we  must  perform  about  16. 1015  multiplications  and 
the  same  number  of  divisions  (which  are  at  least  as  time-consuming).  These 
operation  counts  are  reduced  in  the  author's  program  by  evaluating  and  storing 
some  of  the  quantities  (kgj /(kg j  "  k^))  beforehand,  but  only  by  a  factor  of 
two  (no  greater  saving  than  this  is  possible,  because  the  jC^  term  of  the 
sum  defining  a  u.  , ,  depends  on  all  three  indices,  i,  i'  and  j).  Consequently 
we  could  only  reduce  the  running  time  of  the  program  by  a  small  fraction  even 
if  we  reduced  to  negligible  proportions  Che  number  of  arithmetic  operations 
required  to  solve  the  linear  equations,  so  no  special  algorithms  were  used. 

Finally,  we  consider  the  precision  to  which  the  calculations  must  be 
carried  out.  The  description  in  Somlo’s  paper  [6]  and  his  program  make  it 
clear  that  the  results  in  [  6]  were  calculated  using  about  17  decimal  places 
("double-precision").  The  author's  experience,  however,  is  that  11  decimal 
places  (corresponding  to  the  storing  of  a  real  variable  in  48  bits)  is 
sufficient;  this  is  confirmed  by  the  work  reported  in  [  9) .  Single-precision 
arithmetic  gave  this  number  of  places  on  the  ICL-1906S  computer  used  by  the 
author,  which  saved  substantially  on  time  and  storage. 

20  ERRORS  AND  TOLERANCES  IN  PRACTICAL  APPLICATION 

The  work  described  above  was  undertaken  with  the  intention  of  using 
calculable  coaxial-line  discontinuities  to  construct  calculable  reflection  and 
transmission  standards.  We  therefore  discuss  briefly  the  errors  likely  to 
arise  in  such  applications,  considering  (to  fix  ideas)  the  device  shown  in 
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Fig  4,  which  contains  a  section  of  oversize  inrcr  conductor  and  has  normal- 
size  inner  on  either  side  of  this  section.  Devices  of  this  kind  are  already 
in  use  as  semi-calculable  standards,  both  in  the  Microwave  Standards  Division 
of  RSRE  and  at  other  standards  laborato'ies,  and  it  would  be  useful  to  make 
them  fully  calculable  by  doing  away  wit i  the  need  for  approximate  theory  of 
the  discontinuity  effects. 

Assuming  first  that  the  dimensions  of  the  device  are  exactly  known,  we 
remark  that,  since  the  errors  in  the  capacitances  at  each  of  its  two  steps 
are  a  few  parts  in  10^,  and  since  the  capacitances  are  responsible  for  about 
10%  of  the  total  reflection  (this  was  established  by  calculating  the  reflection 
with  and  without  allowance  for  capacitance),  the  uncertainty  in  reflection 
coefficient  due  to  direct  error  in  the  capacitances  is  a  few  parts  in  10"*. 
Analysis  over  a  wide  frequency  band  is  performed  by  interpolation  between 
capacitance  values  calculated  at  "spot"  frequencies,  but  as  the  frequency 
variation  of  capacitance  is  very  small  (1-2%  between  dc  md  18  GHz  in  standard 
7  mm  line),  the  error  due  to  inaccurate  interpolation  will  be  comparable 
with  the  errors  in  the  "spot"  values.  Finally,  the  effect  of  mutual  interference 
between  the  steps  may  be  assessed  as  negligible.  From  [  1] ,  [2]  and  [  13]  it 
is  known  that  the  presence  of  a  nearby  reflector  at  a  distance  l  from  a  step 
modifies  the  reflection  of  an  evanescent  mode  of  decay  constant  y  by  a  factor 
varying  between  coth  (]y£)  and  tanh  (]y£)  depending  on  t  ie  phase  of  the  TEM 
excitation.  For  small  reflection  these  quantities  are  b >th  nearly  1,  and 

*  •  •  — V  .  , 

their  difference  is  4e  .  If  this  difference  is  negligible  for  even  the 
smallest  significant  y  (which  may  be  obtained  from  the  uoper  critical  frequency 
discussed  earlier)  then  interference  can  be  ignored.  As  an  example,  the  upper 
critical  frequency  for  7  mm  line  has  already  been  determined  (see  section  6, 
"Generalised  modes  and  critical  frequencies")  to  be  75.1  GHz,  so  k'  of  eqn 
(2d)  is  (2ff/c)  x  75.1  -  1.5  mm  Up  to  18  GHz  the  difference  between 

k'  and  (k1  -  k  )  is  fairly  small,  so  we  may  take  y  ^  k'  ^  1.5  mm  .  It  is 
then  clear  that  £  need  only  be  a  few  millimetres  to  make  4e  an  insignificant 
quantity;  and  the  total  error  due  to  im;  erfections  in  calculating  the 

capacitances  will  simply  be  the  sum  of  the  interpolation  and  direct  calculation 

.  .  4 

errors  discussed  above,  about  1  m  10  . 

The  situation  is  quite  different  when  we  come  to  consider  the  effect 
of  imperfectly-known  dimension;.  This  has  been  examined  directly;  four 
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devices  of  the  type  shown  in  Fig  4,  based  on  standard  7  mm  line  and  designed 
to  have  convenient  maximum  VSWRs,  were  analysed,  and  the  calculations  were  then 
repeated  with  all  the  radii  altered  (in  the  directions  of  greatest  effect)  by 
0.001  mm,  which  is  the  error  of  th  •  measurements  of  radius  as  quoted  by  the 
RSRE  metrology  section.  The  general  features  of  the  frequency  variation  of 
the  reflection  coefficient  of  the  device  in  Fig  4  are  shown  in  Fig  5,  and  it 
is  found  that  the  predicted  values  of  the  reflection  coefficient's  maxima  may 
change  by  as  much  as  0.002  for  typical  sizes  of  oversized  inner  conductor. 

These  changes  are  small  enough  to  give  acceptably  small  uncertainties  on  the 
predicted  values  (less  than  0.003  total  uncertainty),  but  when  they  are  com¬ 
pared  with  the  uncertainties  due  to  intrinsic  error  in  the  calculated 
capacitances,  it  is  clear  that  the  precision  of  the  predictions  is  limited  by 
the  mechanical  tolerances;  the  error  of  calculation  is  at  least  an  order  of 
magnitude  smaller.  A  similar,  but  even  stronger,  conclusion  is  reached  in  the 
analysis  of  devices  like  the  one  in  Fig  6,  whose  general  behaviour  is  shown 
in  Fig  7.  The  calculations  for  these  were  performed  using  a  novel  method  of 
accurately  analysing  multi-port  networks  and  cascaded  two-ports,  which  will  be 
described  in  a  future  paper.  Compar  sons  of  measurement  and  prediction  for 
some  of  the  devices  represented  in  F  gs  4  and  6,  with  a  fuller  error  analysis, 
will  also  be  presented  at  an  early  future  date. 

21  CONCLUSIONS 

It  has  been  shown  that  equivalent  circuits  for  problems  involving  junctions 
of  two  generalised  cylindrical  waveguides  can  be  found,  and  a  method  of  cal¬ 
culating  their  component  values  has  been  fully  described  for  the  important 
special  case  of  a  coaxial  line  with  a  radial  step  on  its  inner  conductor 
(including  the  degenerate  form  of  this  problem  in  which  the  inner  conductor  is 
truncated).  The  method  has  been  programmed  in  Algol-68R  on  an  ICL-1906S 

computer,  for  both  the  normal  and  degenerate  forms  of  the  problem,  and  gives 

4 

results  agreeirg  to  a  few  parts  in  1C  with  independent  published  work.  The 
results  are  being  applied  to  the  developmeit  of  calculable  reflection  and 
transmission  standards  for  microwave  aeasui  ements  using  coaxial  lines;  in  this 
application  it  has  b  ;en  shown  that  th ;  accuracy  of  prediction  of  performance 
is  limited  in  practi  :e  only  by  mechanical  tolerances. 
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